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In the 1980s an important goal of the emergent held of fractals was to determine the 
relationships between their physical and geometrical properties. The fractal-Einstein 
and Alexander-Orbach laws, which interrelate electrical, diffusive and fractal prop- 
erties, are two key theories of this type. Here we settle a long standing controversy 
about their exactness by showing that the properties of a class of fractal trees violate 
both laws. A new formula is derived which unifies the two classical results by proving 
that if one holds, then so must the other, and resolves a puzzling discrepancy in the 
properties of Eden trees and diffusion limited aggregates. The failure of the classical 
laws is attributed to anisotropic exploration of the network by a random walker. The 
occurrence of this newly revealed behaviour means that numerous theories, such as 
recent hrst passage time results, are restricted to a narrower range of networks than 
previously thought. 



Consider two simple experiments performed on an arbitrary network of sites linked by 
bonds of identical length. First, release a vast number of random walkers at the centre and 
measure their average distance from the origin (r). Then replace the bonds by resistors, 
apply a unit voltage at a point, and earth all of the sites on a sphere centred at that point. 
The current gives the electrical resistance p of the network. Although simple to conceive, 
the resistance p and distance (r) are fundamental properties: in addition to quantifying 



mass and electronic transport in materials jl, |2j, |3|, |4], |5|, |6|, |7|, |8|,l9|, |10|, [llj, they can be linked 



to a variety of problems such as oil recovery in porous rocks 2] , chemical reaction rates 
and cellular processes [4j, |5| . These elementary properties are also connected to first passage 
times on networks, which have been used to model processes as diverse as viral infections 

n 

and animal foraging strategies [3J. 

A ubiquitous feature of fractals is that their properties follow power laws with non-integer 
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exponents [l|,ll2j]. For example, the mass within a radius r scales as M(r) ~ r d f ; the distance 
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FIG. 1: A model fractal tree, (a) The asymmetric base unit; (b) the second generation of the 
network obtained by attaching multiple copies of the rescaled base unit to the end points (1) and 
(2) of the base unit; (c) an equivalent resistor network used to calculate p* . 

travelled by a random walker scales with time as (r) ~ t 1 ^; the probability that a random 
walker is at its origin scales as c(0, t) ~ t~ d l 2 ; and the electrical resistivity between two 
points scales as p(r) ~ r^. The exponent dj is the fractal dimension, d w is known as the 
random walk dimension, d the spectral dimension, and ( the resistivity exponent. 

The interrelationships between these exponents are the structure-property correlations for 
fractal networks. The Alexander and Orbach[ll] law states that d = 2df/d w and Rammal 
and Toulouse 2J predicted that £ = df(2 — d)/d. Combining both results gives the fractal 
Einstein law lj ( = d w — df, so-called because it can also be derived from a result due 
to Einstein[l|. Although there is preponderance of evidence in their favour, the exactness 
of both formulae is controversial [8|, |9j: computations on two important fractals appear to 
violate the laws. Jacobs et al[3] have used simulation to show that three-dimensional DLA 
has d ~ 1.35 , whereas 2df/d w pa 1.55. Furthermore, their results predict ( = d w — df pa 0.71, 
which disagrees with the estimated value i], [jjj] £ pa 1. The explanation of these exceptions 
is a long-standing challenge in the field [9]. 

To investigate this problem we have studied the properties of a class of fractal trees 
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18| . an example of which is shown in Figure 1. The network is made by taking a 
base unit, doubling its size, and attaching (i = 1, 2) copies of the re-scaled unit to each of 
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the two end points of the base. Continuing the process indefinitely gives an infinite network 
with df = ln(2(M! + u 2 ))/ln(2) and d w = 2. To find the resistance p* between the origin 
and infinity 15], we represent the infinite network by three resistors: the stem of length 1 
and two branches of resistivity p\ and p 2 . Kirchoff's laws states that p\ = 1 + 2p*/u± and 
P2 = 2 + 2p*/u2 , where 2p* is the resistance of the infinite branches connected to the two 
end points of the base unit. Because the shape of each branch is identical with the original 
infinite network (but each element is twice as long), their resistance is 2p*. Kirchoffs law 
for the three- resistor circuit gives p* = 1 + 1/(1/ p\ + 1/ P2) ? which is quadratic in p*. 
Renormalisation methods can be used to derive (Haynes & Roberts, in preparation) the 
spectral dimension and resistivity exponent. For (ui, u 2 ) 7^ (1,1), the quadratic has a finite 
positive root, and it can be shown that 

I u\u 2 + 4(p* 2 + 3p*u 2 + uf)Ui + 4p* 2 uj I ' w w 

For the case (wi, W2) = (1,1)) the quadratic has no positive solutions (implying p* is infinite), 
and we can show d — In (2 (ui + U2)) / ln(2). For all cases, the resistivity exponent is ( = 
2-d. 

The properties of the network violate the Alexander-Orbach and fractal Einstein rela- 
tionships unless Hi = «2 = 1 or «i = 2u 2 . For example, u\ = 1 and u 2 = 20 gives d ~ 4.37, 
which does not equal 2dj/d w = 5.39. Similarly, ( = 2 — d = —2.37 disagrees with the 
prediction ( = d w — df = —3.34 . Standard computations (Appx. [A]) were used to verify d, 
d w , df, and (. 

In order to derive a new relationship between the electrostatic and diffusive properties of 
a network, consider the concentration field generated by the release of a random walker at 
the origin at every time step. This concentration is exactly given by C{r, i) = J °° c(r, r)dr, 
where c(r, t) is the probability of finding a random walker at r, after time t, if a single walker 
is released at the origin at t = 0. To link the dynamic and static problems the integration 
is terminated at T = (R/b) dw , where b is a number of order one. As R = bT l l dw is a 
typical distance reached by the initial walker after time T, only a very small proportion 
of the T + 1 walkers released will exceed this radius; hence C(r,T) m for r > R. In 
the central region the spatial concentration profile C (r, T) is assumed to have equilibrated, 
and therefore satisfies the potential equation. The boundary conditions correspond to the 
potential on a finite network grounded at radius R due to the supply of unit current at the 
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TABLE I: Equation (2) provides a significantly better estimate of £ than the fractal Einstein 
relationship for Eden trees and three dimensional DLA clusters. The table is adapted from Ref. [j| 



using data from Refs 



. 0, 0, 
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Fractal 


d 


d w 


df 


C d m i n 


C = d w -d f 


Eq. (2) 


DLA 2D 


1.20±0.05[13] 


2.64 ± 0.05 [12] 


1.70±0.02[13] 


1.00 ± 0.02 [10] 


0.94 ± 0.07 


1.05 ± 0.09 


DLA 3D 


1.35±0.05[13] 


3.19 ± 0.08[13J 


2.48±0.02[13] 


1.02 ± 0.03 


0.71 ±0.10 


1.04 ±0.10 


Eden tree 2D 


1.22±0.02[20] 


2.82±0.06[21] 


2 


1.22 ±0.02(20] 


0.82 ± 0.07 


1.10 ±0.05 


Eden tree 3D 


1.32±0.02[20] 


3.85 ±0.15[21] 


3 


1.32 ±0.02(20] 


0.85 ± 0.07 


1.31 ±0.09 



origin (I = 1). Now the resistance is simply p(R) ~ C(0, (R/b) dw )/I, which implies 

/•(R/fe) d » 

p(i?) « / c(0,t)tft ~ < 



\og(R) d = 2 

Q R d w (2-d)/2 J > 2 



where Q is a constant. This exactly matches the known scaling behaviour of the resistance 
if 

( = d w (2-d)/2. (2) 



Note that the spectral dimension 19], and hence ( , are site independent, even though p* 
can vary from site to site if d > 2 . Equation (2) is exact for the fractal trees depicted in 
Fig. 1 as d w = 2. Computational checks on the result are provided in Appx. lAl 



Table I shows available simulated data [9J for the properties of Eden trees 20L 21 1 and 

nn 

DLA clusters [lOj, [13J. The resistance of loopless fractals is proportional to the length of the 



shortest path £ between two sites which scales [9[ as £ ~ tf m,n ( so ( = d min ). Eq. (2) 
is seen to provide a good estimate of £ for DLA and Eden trees in three dimensions (thus 
contrasting with the fractal Einstein relation). In two dimensions, Eq. (2) is superior to 
the fractal Einstein relation for Eden trees, whereas for DLA both Eq. (2) and the fractal 
Einstein relation have a similar level of accuracy and are consistent with £ = 1. Data for 
Eden trees were obtained for relatively small clusters, and it would be useful to reconsider 
the calculations. 

The fractal Einstein formula has been rigorously proven [3] using certain assumptions 
about the geometry and a technical "smoothness" criterion on the electrostatic potential. In 



1995, Teles proved 



231 ] that ii ( = d w — df then d = 2df / d w for loopless "smooth" networks 
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with ( > 0; however, the precise connection for general networks has not been established. 
The result in Eq. (2) provides this link: it shows that d = 2df/d w if, and only if, ( = d w — df, 
irrespective of the presence of loops or the sign of £. 

A technical explanation of the assumptions underlying the Alexander-Orbach and fractal 
Einstein laws, and why they do not hold for the fractal tree, is given in Appx. [B]). In 
summary, both results implicitly assume that all sites at a distance r from the origin are 
approximately explored uniformly by a random walker. Although the network may be 
spatially anisotropic, the probability fields, and hence electrostatic fields, are isotropic on 
the network. For the fractal tree which violates both laws, the random walker probability 
(and hence potential and current) fields are anisotropic on the network. This implies that 
a random walker shows a preference for certain directions. Eq. (2) does not require an 
isotropy condition because the electrostatic and probability fields (and hence (, d w and d) 
are similarly affected by anisotropy (Appx. 03]). Figure 2 depicts the current and potential 
distribution on a DLA cluster. Our results indicate that the failure of the Alexander-Orbach 
and fractal Einstein laws are linked to the dramatic non-uniformities in the potential field 
and current distribution of the network. 

In summary, we derive a third fundamental law interrelating the properties of networks. 
This establishes a simple and direct connection between the Alexander-Orbach and fractal 
Einstein laws; if one holds, then so does the other. Because the derivation does not require 
the network to have a fractal dimension, we propose that the relationship holds for inhomo- 
geneous networks 24j (i.e., networks which have no df). Two examples are provided in the 
Appx. El). Note that many results (e.g. refs. [si, for fractal networks are derived using 
the fractal Einstein or Alexander-Orbach formula. This limits their application to networks, 
such as percolation clusters or Sierpinski graphs, on which these relationships are known to 
hold. In addition to advancing our understanding of diffusion on fractals, we believe Eq. (2) 
will find direct application in the fields of science that rely on fractal network models. 
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FIG. 2: Electrostatic fields on a diffusion limited aggregation cluster. The gray branches carry 
current and the potential of the non-current carrying branches is given by their colour. The red 
branches, which emanate near the origin, have the highest potential. Outer branches with potential 
less than 10% of the maximum have been eliminated for clarity. The image was produced using 
specialised software 26]. 
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FIG. 3: The spectral and random walk dimensions. The probability that a walker revisits its origin 
at time t (left hand axis); we calculate d = 4.37 (Best Fit: 4.36). On the right hand axis we plot 
the characteristic distance travelled by the walker \Jr 2 (t) ~ t l l dw against i 1 / 2 confirming d w = 2. 

APPENDIX A: COMPUTATIONAL VERIFICATION OF OUR PRINCIPAL 

RESULTS 

Computations of d and d w are depicted in Figure 1. The results are in excellent agreement 
with the analytic results. To verify Equation (2) p(R) was computed directly for the cases 
(ui,U2) = (1,20) and {u\,U2) = (1,3). Fits to the plots (Fig. 2) of p(R) are in very good 
agreement with the analytic result ( = 2 — d. The scaling relationship is also confirmed on a 
network with ( > obtained by taking (ui, u 2 ) = (1, 1) and quadrupling the branch lengths 
at each iteration (Haynes and Roberts, in preparation). In addition Eq.(2) is analytically 
confirmed (Haynes and Roberts, in preparation) for a non-trivial generalisation of the tree 
obtained by replacing the four elemental branches (of length 2 n ) of the nth generation of 
the network by a deterministic tree|9] of iteration n (using a T shaped generator). In this 
case, C is unchanged because the dangling branches of the T-fractal do not contribute to the 
resistance, but both d w = ln(6)/ln(2) and d ^ 2df/d w differ from the original tree. 
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FIG. 4: Comparison of electrostatic and diffusive properties. The computed resistivity p(R) of three 
networks (symbols) alongside the concentration (7(0, (R/b) dw ) (dashed lines) at the origin of the 
constant source diffusion problem. The value of b (chosen by eye) shifts the lines horizontally. The 
solid lines represent lines of best fit to p(R) : (o) The network with (ui,^) = (1,20), £ = —2.37, 
(Best Fit:-2.28 ), b = 0.9 ; (A) The network with («i,u 2 ) = (1,3), C = -0.941 (Best Fit:-0.92), 
b = 1.1 ; (□) A network with C > discussed in the text, ( = 0.5 (Best Fit:0.527 ), b = 1.1 . 

APPENDIX B: ISOTROPY AND ANISOTROPY ON THE NETWORK 



In the main text it is demonstrated that the Alexander- Orbach (AO) and fractal Einstein 
(FE) laws do not apply to certain fractal networks. A new accurate formula for the resistivity 
exponent is derived. Both findings can be explained heuristically by considering the role of 
anisotropy on the network. 

The probability that a walker released at r = at t = will be at the point r on a 
network after time t will depend on the direction of r as well as its magnitude r = |r|. This 
probability is denoted as c a (r, t), where the subscript a (anisotropic) differentiates it from 
the function c(r,t) used in the main text. The two functions can be related by 

c(r,t) = -i- [ c a (r,t)dS, (Bl) 
S{r) J s 

where S(r) ~ r d ^ x is the area (mass) of the fractal at radius r. c(r,t) can be regarded as 
a network-spherical average, because the average on the shell is only taken over the regions 
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occupied by the network. Equivalently, it can be called 25[ the average probability per site. 

After time t, a walker released from the origin will on average have explored a region of 
radius R ~ t 1 ' dw . As there are V(r) ~ R df sites within that radius, f v , R \ c (r, t)dV = P 
where P is the probability that the walker is in the central region. If this region is explored 
approximately uniformly [9], then c (r, t)/c a (0, t) is a slowly varying function for r < R. The 
long time behaviour of the probability is found by setting c (r, t) ~ c o (0, t) which gives 
c a (0,t) J V(R) dV w P or c(0,t) « P/V(R{t)) ~ This provides the rationale behind 

the AO law d = 2df/d w . The derivation assumes that the volume V(r(t)) is approximately 
uniformly explored for |r| < i?. In particular, this requires that c a (r,t) ~ c(r, t), i.e., 
the concentration field is approximately isotropic on the network (or network-spherically 
symmetric). If this is not true, the volume explored by the walker will generally not be 
V(R(t)). Data shown in Fig. [5] confirms that c a (r,t) « c(r,t) for the fractal tree with 
(ui, u 2 ) = (1,2). This tree's properties follow the AO law exactly. In contrast c a (r, t) is seen 
to be strongly anisotropic for the fractal tree with (^1,1*2) = (20, 1). The breakdown of the 
AO law is attributed to the non-uniform exploration of the network. For these classes of 
fractal trees, the validity of the AO law is seen to depend on the probability density being 
isotropic on the network. 

A similar requirement of uniformity is implicitly assumed in the derivation of the FE law. 
This is clearly seen in an examination of the total current flow / through a shell of thickness 
Ar. By definition I = ne x dS x Ar/(At) where n is the charge carrier density, e is the 
carrier charge, dS is an element of area, and At is the time it takes a charge to cross the 
shell. Now the time scale for diffusing a distance r is t ~ r dw , so At ~ r*"~ :L Ar. Summing 
over the total area of the shell gives 

/ „ Ar / ^ « Ar4- ~ r d ^ dw , 
Js At At 

and therefore p = V/I ~ r d f~ dw which proves the FE law. Although the argument assumes 
At is independent of direction, this is only strictly true if c a (r, t) is uniform over the shell. 
Therefore, if diffusion exhibits preferential directions on the network (as it does for the tree 
with (1*1,^2) = (20, 1)), the FE law will be invalid. 

The fractal tree provides a concrete example of the qualitative balance arguments ex- 
pressed above. As noted in the main text, the FE law is obeyed for the case U\ = 2u 2 - 
Rearranging the expressions for pi and p 2 gives p 2 u 2 — p\U\ = 2u 2 — u\\ hence the condition 
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FIG. 5: The spherically averaged functions c(r,t) for the fractal tree with (uijUz) = (1,2) (a) and 
(u\,U2) = (20,1) (□). The lines indicate the highest and lowest values of c a (r,t) at each r. The 
significant anistropy in the latter model is connected to the breakdown of the AO and FE laws. 

u\ = 2«2 implies p<il P\ = h/h = U\/u%, where Ii (i = 1,2) are the currents on each branch 
respectively. As the ratio of the masses of the branches is U\jui , it is seen that conventional 
scaling holds because the mass and current on different branches extending from a node 
are balanced. However, for the case (111,1x2) = (1,20), a significant mass-current inbalance 
occurs; although there is twenty times more mass in branch 1 than branch 2, the current is 
only about /1//2 ~ 3 times greater. Although this discussion provides a useful picture of 



why the FE law fails for fractal trees studied here, note that Teles 22| has provided general 
conditions in terms of the potential field. 

This simple numerical example also illustrates the breakdown of the AO law in terms of 
a mass-probability imbalance. It is known that the ratio of currents I1/I2 = P1/P2, where 
Pi is the probability that a walker on branch % never returns to its origin (i.e., it escapes). 
Consider a walker at the first junction in the network. There is 20 times more mass in 
branch 1 than branch 2, but the probability of escaping along branch 1, and never returning 
to the junction, is only about 3 times greater than that of escaping along branch 2. The 
greater-than-expected return of walkers from branch 1 increases c(0, t), and qualitatively 
explains why d < 2df/d w . 
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APPENDIX C: INHOMOGENEOUS NETWORKS 



Some classes 24] of important networks do not have a fractal dimension df (called inho- 
mogeneous networks). As the derivation of Eq. (2) does not directly involve the mass, we 
hypothesise that it can be used to predict the properties of these networks. Technically, d w 
only exists for fractal networks, but an analogous exponent a is defined by < r 2 >~ t a , 
whereby Eq. (2) becomes ( = (2 — d)/a. Consider the fractal tree, which becomes inhomo- 
geneous as w 2 — > 00 ( dtf ~~ * °° ) • The argument used to show ( = 2 — d is not altered in this 



limit, so Eq. (2) continues to hold. A second example is provided by the comb lattice 1] with 
an infinite spine and infinite teeth. This inhomogeneous network has d = 3/2 and a = 1 
(since the teeth are one- dimensional). The exponent ( can be easily found if the teeth are 
folded against the spine (This will not affect d and a). If the lattice is earthed along a line 
perpendicular to the spine (at node 0), the resistance between the nth node to the left of 
the line and the line is p n = l/(l/n + l/p n _i) with pi = 1/2. For large n p n ~ p n _i, and 
the solution of the quadratic is p n = \fn so ( = 1/2 which is consistent with ( = (2 — d)/a. 



APPENDIX D: AN EXPLANATION OF EQ. (2) 

It is instructive to consider the relationship between the mean-squared distance and 
potential in terms of their connection with the anisotropic probability density c a (r, t). There 
are numerous ways of defining resistance on a network. The point-to-shell resistance is 
defined by earthing all sites a distance R from an origin, and applying a potential at that 
origin. The resultant potential field is denoted as (f)(r;R), hence the resistance is p(R) = 
0(0; R)/I where / is the current. By the argument given in the main text, the potential at 
a site r is 

r(R/b) d w 

a (r;i2)« / c a (r,t)dt. (Dl) 

Taking the network average 



o 



s 



1 r I r(R/b) dw r 

- / M^R)dS &g / c a (r,t)dSdt (D2) 



and applying the scaling arguments gives Eq. (2) of the main text 



C = y(2-rf). (D3) 
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The mean squared distance is directly related to the probability density by 

r poo r 

(r 2 (t))= / c (r, t)r 2 dV = r 2 c a (r,t)dSdr. (D4) 
Jv Jo Js 

We believe that Eqn. (1D3I) can hold for structures which do not adhere to the AO and 

FE laws because it does not depend explicitly on the distribution of mass M(r). Consider 

the three exponents d, d w (or a) and £. The probability c(0, t) ~ t~ d l 2 does not depend on 

the distribution of M(r) because the network can be distorted arbitrarily (if the inter-site 

distance is preserved) without changing c(0, t) and hence d. The exponents d w (or a) and £ 

may be affected by a spatial distortion of the network, but as both variables are ultimately 

linked to spherical averages of c a (r, t) by Eqn. (ID2I) or flD4j) . they will change in a consistent 

manner. In terms of the comb example given above; variation of the angle between the teeth 

and spine does not alter c(0, t) and will affect c a (r, t) and a (r, t) similarly. For DLA clusters, 

and the fractal tree that violates the AO and FE laws, the fact that massive portions of 

the structure are not in balance with the distribution of potential or probability does not 

disturb the relationship between d, d w and (. 



[1] Havlin, S. Sz ben Avraham, D. Diffusion in disordered media. Adv. Phys. 51, 187 (2002). 
[2] Sahimi, M. Flow phenomena in rocks: from continuum models to fractals, percolation, cellular 

automata, and simulated annealing. Rev. Mod. Phys. 65, 1393-1534 (1993). 
[3] Condamin, S., Benichou, O., Tejedor, V., Voituriez, R. & Klafter, J. First-passage times in 

complex scale-invariant media. Nature 450, 77 (2007). 
[4] Gallos, L. K., Song, C, Havlin, S. &: Makse, H. A. Scaling theory of transport in complex 

biological networks. Proc. Natl. Acad. Sci 104, 7746-7751 (2007). 
[5] Song, C. M., Havlin, S. & Makse, H. A. Self-similarity of complex networks. Nature 433, 392 

(2005). 

[6] Orbach, R. Dynamics of fractal networks. Science 231, 814-819 (1986). 

[7] Rammal, R. &: Toulouse, G. Random-walks on fractal structures and percolation clusters. J. 

Phys. Lett. (Paris) 44, L13 (1983). 
[8] Hughes, B. D. Random walks and random environments. Volume 2 (Clareton Press, Oxford, 

1996). 



13 



[9] ben Avraham, D. & Havlin, S. Diffusion and reactions in fractals and disordered systems 
(Cambridge Univ. Press, Cambridge, uK, 2000). 
[10] Meakin, P., Majid, L, Havlin, S., Stanley, H. E. k Witten, T. E. Topological properties 
of diffusion limited aggregation and cluster-cluster aggregation. J. Phys. A 17, L975-L981 
(1984). 

[11] Alexander, S. & Orbach, R. Density of states on fractals: fractons. J. Phys. (Paris) Lett. 19, 
L625 (1982). 

[12] Mandelbrot, B. B. The fractal geometry of nature (Freeman, San Francisco, 1982). 

[13] Jacobs, D. J., Mukherjee, S. & Nakanishi, H. Diffusion on a DLA cluster in two and three 

dimensions. J. Phys. A: Math. Gen. 27, 4341-4350 (1994). 
[14] Witten, T. A. & Sander, L. M. Diffusion-limited aggregation, a kinetic critical phenomenon. 

Phys. Rev. Lett. 47, 1400-1403 (1981). 
[15] Doyle, P. G. & Snell, J. L. Random walks and electric networks. (Math. Assoc. of America, 

Washington, DC, USA, 1984). 
[16] Burioni, R. & Cassi, D. Spectral dimension of fractal trees. Phys. Rev. E 51, 2865 (1995). 
[17] Kron, B. & Teufl, E. Asymptotics of the transition probabilities of the simple random walk 

on self-similar graphs. Trans. Amer. Math. Soc. 356, 393 (2004). 
[18] Haynes, C. P. & Roberts, A. P. Spectral dimension of fractal trees. Phys. Rev. E (accepted) 

(2008). 

[19] Hattori, K., Hattori, T. & Watanabe, H. Gaussian field theories on general networks and the 

spectral dimensions. Progr. Theoret. Phys. (Suppl.) 92, 108 (1987). 
[20] Nakanishi, H. & Herrmann., H. J. Diffusion and spectral dimension on Eden tree. J. Phys. 

A: Math. Gen. 26, 4513-4519 (1993). 
[21] Reis, F. D. A. A. Scaling for random walks on Eden trees. Phys. Rev. E. 92, 3079-3081 

(1996). 

[22] Teles, A. Random walks on graphs, electric networks and fractals. Probab. Th. Rel. Fields 
82, 435 (1989). 

[23] Teles, A. Spectra of graphs and fractal dimensions II. J. Theoret. Probab. 8, 77 (1995). 
[24] Cassi, D. & Regina, S. Random walks on bundled structures. Phys. Rev. Lett. 76, 2914 (1996). 
[25] O'Shaughnessy, B. & Procaccia, I. Analytical solutions for diffusion on fractal objects. Phys. 
Rev. Lett. 54, 455 (1985). 



14 

[26] Bourke, P. B. Constrained diffusion limited aggregation in 3 dimensions. Comp. Graphics. 
30, 640-649 (2006). 



Acknowledgements Paul Bourke, University of Western Australia, produced the 
DLA visualisations from the authors' data. 



